Genetic heterogeneity in hereditary hearing loss: Potential role of kinociliary protein TOGARAM2.

Hearing loss (HL) is a heterogenous trait with pathogenic variants in more than 200 genes that have been discovered in studies involving small and large HL families. Over one-third of families with hereditary HL remain etiologically undiagnosed after screening for mutations in the recognized genes. Genetic heterogeneity complicates the analysis in multiplex families where variants in more than one gene can be causal in different individuals even in the same sibship. We employed exome or genome sequencing in at least two affected individuals with congenital or prelingual-onset, severe to profound, non-syndromic, bilateral sensorineural HL from four multiplex families. Bioinformatic analysis was performed to identify variants in known and candidate deafness genes. Our results show that in these four families, variants in a single HL gene do not explain HL in all affected family members, and variants in another known or candidate HL gene were detected to clarify HL in the entire family. We also present a variant in TOGARAM2 as a potential cause underlying autosomal recessive non-syndromic HL by showing its presence in a family with HL, its expression in the cochlea and the localization of the protein to cochlear hair cells. Conclusively, analyzing all affected family members separately can serve as a good source for the identification of variants in known and novel candidate genes for HL.

Identification of novel MYH14 variants in families with autosomal dominant sensorineural hearing loss.

Autosomal dominant sensorineural hearing loss (ADSNHL) is a genetically heterogeneous disorder caused by pathogenic variants in various genes, including MYH14. However, the interpretation of pathogenicity for MYH14 variants remains a challenge due to incomplete penetrance and the lack of functional studies and large families. In this study, we performed exome sequencing in six unrelated families with ADSNHL and identified five MYH14 variants, including three novel variants. Two of the novel variants, c.571G > C (p.Asp191His) and c.571G > A (p.Asp191Asn), were classified as likely pathogenic using ACMG and Hearing Loss Expert panel guidelines. In silico modeling demonstrated that these variants, along with p.Gly1794Arg, can alter protein stability and interactions among neighboring molecules. Our findings suggest that MYH14 causative variants may be more contributory and emphasize the importance of considering this gene in patients with nonsyndromic mainly post-lingual severe form of hearing loss. However, further functional studies are needed to confirm the pathogenicity of these variants.

Genome sequencing identifies coding and non-coding variants for non-syndromic hearing loss.

Hearing loss (HL) is a common heterogeneous trait that involves variants in more than 200 genes. In this study, we utilized exome (ES) and genome sequencing (GS) to effectively identify the genetic cause of presumably non-syndromic HL in 322 families from South and West Asia and Latin America. Biallelic GJB2 variants were identified in 58 probands at the time of enrollment these probands were excluded. In addition, upon review of phenotypic findings, 38/322 probands were excluded based on syndromic findings at the time of ascertainment and no further evaluation was performed on those samples. We performed ES as a primary diagnostic tool on one or two affected individuals from 212/226 families. Via ES we detected a total of 78 variants in 30 genes and showed their co-segregation with HL in 71 affected families. Most of the variants were frameshift or missense and affected individuals were either homozygous or compound heterozygous in their respective families. We employed GS as a primary test on a subset of 14 families and a secondary tool on 22 families which were unsolved by ES. Although the cumulative detection rate of causal variants by ES and GS is 40% (89/226), GS alone has led to a molecular diagnosis in 7 of 14 families as the primary tool and 5 of 22 families as the secondary test. GS successfully identified variants present in deep intronic or complex regions not detectable by ES.

H4C5 missense variant leads to a neurodevelopmental phenotype overlapping with Angelman syndrome.

Recurrent de novo missense variants in H4 histone genes have recently been associated with a novel neurodevelopmental syndrome that is characterized by intellectual disability and developmental delay as well as more variable findings that include short stature, microcephaly, and facial dysmorphisms. A 4-year-old male with autism, developmental delay, microcephaly, and a happy demeanor underwent evaluation through the Undiagnosed Disease Network. He was clinically suspected to have Angelman syndrome; however, molecular testing was negative. Genome sequencing identified the H4 histone gene variant H4C5 NM_003545.4: c.295T>C, p.Tyr99His, which parental testing confirmed to be de novo. The variant met criteria for a likely pathogenic classification and is one of the seven known disease-causing missense variants in H4C5. A comparison of our proband's findings to the initial description of the H4-associated neurodevelopmental syndrome demonstrates that his phenotype closely matches the spectrum of those reported among the 29 affected individuals. As such, this report corroborates the delineation of neurodevelopmental syndrome caused by de novo missense H4 gene variants. Moreover, it suggests that cases of clinically suspected Angelman syndrome without molecular confirmation should undergo exome or genome sequencing, as novel neurodevelopmental syndromes with phenotypes overlapping with Angelman continue to be discovered.

Mutations in MINAR2 encoding membrane integral NOTCH2-associated receptor 2 cause deafness in humans and mice.

Discovery of deafness genes and elucidating their functions have substantially contributed to our understanding of hearing physiology and its pathologies. Here we report on DNA variants in MINAR2, encoding membrane integral NOTCH2-associated receptor 2, in four families underlying autosomal recessive nonsyndromic deafness. Neurologic evaluation of affected individuals at ages ranging from 4 to 80 y old does not show additional abnormalities. MINAR2 is a recently annotated gene with limited functional understanding. We detected three MINAR2 variants, c.144G > A (p.Trp48*), c.412_419delCGGTTTTG (p.Arg138Valfs*10), and c.393G > T, in 13 individuals with congenital- or prelingual-onset severe-to-profound sensorineural hearing loss (HL). The c.393G > T variant is shown to disrupt a splice donor site. We show that Minar2 is expressed in the mouse inner ear, with the protein localizing mainly in the hair cells, spiral ganglia, the spiral limbus, and the stria vascularis. Mice with loss of function of the Minar2 protein (Minar2tm1b/tm1b) present with rapidly progressive sensorineural HL associated with a reduction in outer hair cell stereocilia in the shortest row and degeneration of hair cells at a later age. We conclude that MINAR2 is essential for hearing in humans and mice and its disruption leads to sensorineural HL. Progressive HL observed in mice and in some affected individuals and as well as relative preservation of hair cells provides an opportunity to interfere with HL using genetic therapies.

Genome sequencing reveals novel noncoding variants in PLA2G6 and LMNB1 causing progressive neurologic disease.

Neurodegenerative disorders and leukodystrophies are progressive neurologic conditions that can occur following the disruption of intricately coordinated patterns of gene expression. Exome sequencing has been adopted as an effective diagnostic tool for determining the underlying genetic etiology of Mendelian neurologic disorders, however genome sequencing offer advantages in its ability to identify and characterize copy number, structural, and sequence variants in noncoding regions. Genome sequencing from peripheral leukocytes was performed on two patients with progressive neurologic disease of unknown etiology following negative genetic investigations including exome sequencing. RNA sequencing from peripheral blood was performed to determine gene expression patterns in one of the patients. Potential causative variants were matched to the patients' clinical presentation. The first proband was found to be heterozygous for a likely pathogenic missense variant in PLA2G6 (c.386T>C; p.Leu129Pro) and have an additional deep intronic variant in PLA2G6 (c.2035-926G>A). RNA sequencing indicated this latter variant created a splice acceptor site leading to the incorporation of a pseudo-exon introducing a premature termination codon. The second proband was heterozygous for a 261 kb deletion upstream of LMNB1 that included an enhancer region. Previous reports of copy number variants spanning this region of cis-acting regulatory elements corroborated its pathogenicity. When combined with clinical presentations, these findings led to a definitive diagnosis of autosomal recessive infantile neuroaxonal dystrophy and autosomal dominant adult-onset demyelinating leukodystrophy, respectively. In patients with progressive neurologic disease of unknown etiology, genome sequencing with the addition of RNA analysis where appropriate should be considered for the identification of causative noncoding pathogenic variants.

Genome-Wide Linkage and Association Study of Childhood Gender Nonconformity in Males.

Male sexual orientation is influenced by environmental and complex genetic factors. Childhood gender nonconformity (CGN) is one of the strongest correlates of homosexuality with substantial familiality. We studied brothers in families with two or more homosexual brothers (409 concordant sibling pairs in 384 families, as well as their heterosexual brothers), who self-recalled their CGN. To map loci for CGN, we conducted a genome-wide linkage scan (GWLS) using SNP genotypes. The strongest linkage peaks, each with significant or suggestive two-point LOD scores and multipoint LOD score support, were on chromosomes 5q31 (maximum two-point LOD = 4.45), 6q12 (maximum two-point LOD = 3.64), 7q33 (maximum two-point LOD = 3.09), and 8q24 (maximum two-point LOD = 3.67), with the latter not overlapping with previously reported strongest linkage region for male sexual orientation on pericentromeric chromosome 8. Family-based association analyses were used to identify associated variants in the linkage regions, with a cluster of SNPs (minimum association p = 1.3 × 10-8) found at the 5q31 linkage peak. Genome-wide, clusters of multiple SNPs in the 10-6 to 10-8 p-value range were found at chromosomes 5p13, 5q31, 7q32, 8p22, and 10q23, highlighting glutamate-related genes. This is the first reported GWLS and genome-wide association study on CGN. Further increasing genetic knowledge about CGN and its relationships to male sexual orientation should help advance our understanding of the biology of these associated traits.

Genomic evidence consistent with antagonistic pleiotropy may help explain the evolutionary maintenance of same-sex sexual behaviour in humans.

Human same-sex sexual behaviour (SSB) is heritable, confers no immediately obvious direct reproductive or survival benefit and can divert mating effort from reproductive opportunities. This presents a Darwinian paradox: why has SSB been maintained despite apparent selection against it? We show that genetic effects associated with SSB may, in individuals who only engage in opposite-sex sexual behaviour (OSB individuals), confer a mating advantage. Using results from a recent genome-wide association study of SSB and a new genome-wide association study on number of opposite-sex sexual partners in 358,426 individuals, we show that, among OSB individuals, genetic effects associated with SSB are associated with having more opposite-sex sexual partners. Computer simulations suggest that such a mating advantage for alleles associated with SSB could help explain how it has been evolutionarily maintained. Caveats include the cultural specificity of our UK and US samples, the societal regulation of sexual behaviour in these populations, the difficulty of measuring mating success and the fact that measured variants capture a minority of the total genetic variation in the traits.

Genome-Wide Linkage Study Meta-Analysis of Male Sexual Orientation.

Male sexual orientation is a scientifically and socially important trait shown by family and twin studies to be influenced by environmental and complex genetic factors. Individual genome-wide linkage studies (GWLS) have been conducted, but not jointly analyzed. Two main datasets account for > 90% of the published GWLS concordant sibling pairs on the trait and are jointly analyzed here: MGSOSO (Molecular Genetic Study of Sexual Orientation; 409 concordant sibling pairs in 384 families, Sanders et al. (2015)) and Hamer (155 concordant sibling pairs in 145 families, Mustanski et al. (2005)). We conducted multipoint linkage analyses with Merlin on the datasets separately since they were genotyped differently, integrated genetic marker positions, and combined the resultant LOD (logarithm of the odds) scores at each 1 cM grid position. We continue to find the strongest linkage support at pericentromeric chromosome 8 and chromosome Xq28. We also incorporated the remaining published GWLS dataset (on 55 families) by using meta-analytic approaches on published summary statistics. The meta-analysis has maximized the positional information from GWLS of currently available family resources and can help prioritize findings from genome-wide association studies (GWAS) and other approaches. Although increasing evidence highlights genetic contributions to male sexual orientation, our current understanding of contributory loci is still limited, consistent with the complexity of the trait. Further increasing genetic knowledge about male sexual orientation, especially via large GWAS, should help advance our understanding of the biology of this important trait.

Response to Comment on "Large-scale GWAS reveals insights into the genetic architecture of same-sex sexual behavior".

Hamer et al argue that the variable "ever versus never had a same-sex partner" does not capture the complexity of human sexuality. We agree and said so in our paper. But Hamer et al neglect to mention that we also reported follow-up analyses showing substantial overlap of the genetic influences on our main variable and on more nuanced measures of sexual behavior, attraction, and identity.

A founder noncoding GALT variant interfering with splicing causes galactosemia.

Galactosemia is a rare, treatable hereditary disorder of carbohydrate metabolism. We investigated the etiology of decreased GALT enzyme activity in a cohort of newborns referred by the Florida Newborn Screening Program with no detectable GALT variants in diagnostic molecular tests. Six affected individuals from four families with Guatemalan heritage were included. GALT enzyme activity ranged from 20% to 34% of normal. Clinical findings were unremarkable except for speech delay in two children. Via genome sequencing followed by Sanger confirmation we showed that all affected individuals were homozygous for a deep intronic GALT variant, c.1059+390A>G, which segregated as an autosomal recessive trait in all families. The intronic variant disrupts splicing and leads to a premature termination and is associated with a single haplotype flanking GALT, suggesting a founder effect. In conclusion, we present a deep intronic GALT variant leading to a biochemical variant form of galactosemia. This variant remains undiagnosed until it is specifically targeted in molecular testing.

Long-range cis-regulatory elements controlling GDF6 expression are essential for ear development.

Molecular mechanisms governing the development of the mammalian cochlea, the hearing organ, remain largely unknown. Through genome sequencing in 3 subjects from 2 families with nonsyndromic cochlear aplasia, we identified homozygous 221-kb and 338-kb deletions in a noncoding region on chromosome 8 with an approximately 200-kb overlapping section. Genomic location of the overlapping deleted region started from approximately 350 kb downstream of GDF6, which codes for growth and differentiation factor 6. Otic lineage cells differentiated from induced pluripotent stem cells derived from an affected individual showed reduced expression of GDF6 compared with control cells. Knockout of Gdf6 in a mouse model resulted in cochlear aplasia, closely resembling the human phenotype. We conclude that GDF6 plays a necessary role in early cochlear development controlled by cis-regulatory elements located within an approximately 500-kb region of the genome in humans and that its disruption leads to deafness due to cochlear aplasia.

Spectrum of Genetic Variants Associated with Anterior Segment Dysgenesis in South Florida.

Anterior segment dysgenesis (ASD) comprises a wide spectrum of developmental conditions affecting the cornea, iris, and lens, which may be associated with abnormalities of other organs. To identify disease-causing variants, we performed exome sequencing in 24 South Florida families with ASD. We identified 12 likely causative variants in 10 families (42%), including single nucleotide or small insertion-deletion variants in B3GLCT, BMP4, CYP1B1, FOXC1, FOXE3, GJA1, PXDN, and TP63, and a large copy number variant involving PAX6. Four variants were novel. Each variant was detected only in one family. Likely causative variants were detected in 1 out of 7 black and 9 out of 17 white families. In conclusion, exome sequencing for ASD allows us to identify a wide spectrum of rare DNA variants in South Florida. Further studies will explore missing variants, especially in the black communities.

Large-scale GWAS reveals insights into the genetic architecture of same-sex sexual behavior.

Twin and family studies have shown that same-sex sexual behavior is partly genetically influenced, but previous searches for specific genes involved have been underpowered. We performed a genome-wide association study (GWAS) on 477,522 individuals, revealing five loci significantly associated with same-sex sexual behavior. In aggregate, all tested genetic variants accounted for 8 to 25% of variation in same-sex sexual behavior, only partially overlapped between males and females, and do not allow meaningful prediction of an individual's sexual behavior. Comparing these GWAS results with those for the proportion of same-sex to total number of sexual partners among nonheterosexuals suggests that there is no single continuum from opposite-sex to same-sex sexual behavior. Overall, our findings provide insights into the genetics underlying same-sex sexual behavior and underscore the complexity of sexuality.

Novel variant p.E269K confirms causative role of PLS1 mutations in autosomal dominant hearing loss.

Auditory reception relies on the perception of mechanical stimuli by stereocilia and its conversion to electrochemical signal. Mechanosensory stereocilia are abundant in actin, which provides them with structural conformity necessary for perception of auditory stimuli. Out of three major classes of actin-bundling proteins, plastin 1 encoded by PLS1, is highly expressed in stereocilia and is necessary for their regular maintenance. A missense PLS1 variant associated with autosomal dominant hearing loss (HL) in a small family has recently been reported. Here, we present another PLS1 missense variant, c.805G > A (p.E269K), in a Turkish family with autosomal dominant non-syndromic HL confirming the causative role of PLS1 mutations in HL. We propose that HL due to the p.E269K variant is from the loss of a stable PLS1-ACTB interaction.

A truncating CLDN9 variant is associated with autosomal recessive nonsyndromic hearing loss.

While the importance of tight junctions in hearing is well established, the role of Claudin- 9 (CLDN9), a tight junction protein, in human hearing and deafness has not been explored. Through whole-genome sequencing, we identified a one base pair deletion (c.86delT) in CLDN9 in a consanguineous family from Turkey with autosomal recessive nonsyndromic hearing loss. Three affected members of the family had sensorineural hearing loss (SNHL) ranging from moderate to profound in severity. The variant is predicted to cause a frameshift and produce a truncated protein (p.Leu29ArgfsTer4) in this single-exon gene. It is absent in public databases as well as in over 1000 Turkish individuals, and co-segregates with SNHL in the family. Our in vitro studies demonstrate that the mutant protein does not localize to cell membrane as demonstrated for the wild-type protein. Mice-lacking Cldn9 have been shown to develop SNHL. We conclude that CLDN9 is essential for proper audition in humans and its disruption leads to SNHL in humans.

FOXF2 is required for cochlear development in humans and mice.

Molecular mechanisms governing the development of the human cochlea remain largely unknown. Through genome sequencing, we identified a homozygous FOXF2 variant c.325A>T (p.I109F) in a child with profound sensorineural hearing loss (SNHL) associated with incomplete partition type I anomaly of the cochlea. This variant is not found in public databases or in over 1000 ethnicity-matched control individuals. I109 is a highly conserved residue in the forkhead box (Fox) domain of FOXF2, a member of the Fox protein family of transcription factors that regulate the expression of genes involved in embryogenic development as well as adult life. Our in vitro studies show that the half-life of mutant FOXF2 is reduced compared to that of wild type. Foxf2 is expressed in the cochlea of developing and adult mice. The mouse knockout of Foxf2 shows shortened and malformed cochleae, in addition to altered shape of hair cells with innervation and planar cell polarity defects. Expressions of Eya1 and Pax3, genes essential for cochlear development, are reduced in the cochleae of Foxf2 knockout mice. We conclude that FOXF2 plays a major role in cochlear development and its dysfunction leads to SNHL and developmental anomalies of the cochlea in humans and mice.

Identification of candidate gene FAM183A and novel pathogenic variants in known genes: High genetic heterogeneity for autosomal recessive intellectual disability.

The etiology of intellectual disability (ID) is heterogeneous including a variety of genetic and environmental causes. Historically, most research has not focused on autosomal recessive ID (ARID), which is a significant cause of ID, particularly in areas where parental consanguinity is common. Identification of genetic causes allows for precision diagnosis and improved genetic counseling. We performed whole exome sequencing to 21 Turkish families, seven multiplex and 14 simplex, with nonsyndromic ID. Based on the presence of multiple affected siblings born to unaffected parents and/or shared ancestry, we consider all families as ARID. We revealed the underlying causative variants in seven families in MCPH1 (c.427dupA, p.T143Nfs*5), WDR62 (c.3406C>T, p.R1136*), ASPM (c.5219_5225delGAGGATA, p.R1740Tfs*7), RARS (c.1588A>G, p.T530A), CC2D1A (c.811delG, p.A271Pfs*30), TUSC3 (c.793C>T, p.Q265*) and ZNF335 (c.808C>T, p.R270C and c.3715C>A, p.Q1239K) previously linked with ARID. Besides ARID genes, in one family, affected male siblings were hemizygous for PQBP1 (c.459_462delAGAG, p.R153Sfs*41) and in one family the proband was female and heterozygous for X-chromosomal SLC9A6 (c.1631+1G>A) variant. Each of these variants, except for those in MCPH1 and PQBP1, have not been previously published. Additionally in one family, two affected children were homozygous for the c.377G>A (p.W126*) variant in the FAM183A, a gene not previously associated with ARID. No causative variants were found in the remaining 11 families. A wide variety of variants explain half of families with ARID. FAM183A is a promising novel candidate gene for ARID.